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Abstract 

We report a preliminary measurement of the branching fraction for the decay — > K^K^K^, 
where the mesons are reconstructed through the decay vr^vr". The measurement was 

performed on a sample of 211 x 10^ BB pairs collected by the BABAR detector running on the 
T{4:S) resonance at the PEP-II storage ring. The branching fraction is measured to be 

B{B° K'^K'^K'^) = (6.5 ± 0.8 ± 0.8) x IQ-^, 

where the errors are statistical and systematic, respectively. 
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1 INTRODUCTION 



In this paper we report a measurement of the branching fraction (B) for — > K^K^K^. This 
decay is expected to be penguin dominated; the simplest diagram that can be drawn without 
rescattering is shown in Fig. ^ ^ K^K^K^ is not Cabibbo-suppressed, and so is expected to 
have substantially larger branching fraction than B^ 

In Ref. X the branching fraction B{B^ K'^K'^K^) is related to B{B+ K+K-K+). With 
the assumption of gluonic penguin dominance and the usual assumption B{K^ K^) = B{K^ 
K^) = 0.5, they derive B{B+ K+K-K+) = B{B^ K^K^K^) = m{B^ K^K^K^). 
Using the BABAR and Belle averaged value of K+K-K+) = (29.5 ± 1.8) x 10"^, 

K^K^K^) is expected to be ~ 4 x IQ-^. 




Figure 1: The dominant diagram for the decay B K^K^K^. 



The only previous measurement pi of this branching fraction is from the Belle Collaboration, 
B{B^ K^K^K^) = {4.2tli ± 0.8) x 10"^ based on 78 fb^^ of on-resonance data. The data 
sample used in this paper is 2.7 the size of that data sample, and the B^ K^K^K^ efficiency we 
estimate for our analysis is larger than that in Ref. [3]. 

In this paper we report an inclusive measurement of B^ K^K^K^. In addition to non-resonant 
three-body h — > sqq, (qq = 'ss or dd) gluonic penguin decays, charmless resonant intermediate states 
like B^ foK^ can produce the 3Kg final state. There may also be 6 — > ccs decays that lead to 
the 3K^ final state. The dominant of these is expected to be B^ — > XcO^^^ but its product 

branching fraction ^ is < 0.5 x 10~^ (90% CL), about a factor of ten smaller than that expected 
for B^ K^KgK^. We do not exclude these from this measurement, though we will do a search 
for B^ — > XcoKg as a systematic check. The product branching fraction for B^ — > D^K^ — > 3Kg is 
estimated to be ~ 9 x 10^^ and therefore will be ignored. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
We use 191 fb^^ of data taken at the center-of-mass (CM) energy of the T{AS) resonance (the 
on-resonance data sample). These data correspond to 211 x 10^ BB pairs. 

The BABAR detector is described elsewhere 15). The important parts of the detector for this anal- 
ysis are the charged particle tracking detectors. These consist of five layers of double-sided silicon- 
strip detectors between the beampipe and a 40-layer cylindrical drift chamber, with both axial and 
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small-stereo-angle superlayers. Both detectors are in a 1.5 T solenoidal magnetic field, and pro- 
vide excellent pattern recognition and momentum measurement for reconstruction of vr^vr" 
decays. The electromagnetic calorimeter also contributes to this analysis through the reconstruc- 
tion of neutral particles which are used along with charged tracks not coming from the candidate 
— > K^K^K^ decay to form continuum rejection variables. 

Large samples of Monte Carlo (MC) simulated events are used throughout this analysis, to 
estimate the reconstruction efficiency and to derive parameters used to describe the signal and 
background (BG) shapes in the fit for the signal yield. Except for some parametrized MC samples 
used to validate the fit, all the MC samples were generated with GEANT4 ^ , with the full detector- 
response simulation, and reconstructed with the same programs used for data reconstruction. 

3 CANDIDATE SELECTION 

3.1 A'° RECONSTRUCTION 

All Kg candidates used in this analysis are reconstructed through the decay — > tt+tt". Every pair 
of oppositely-charged tracks that pass a very loose mass selection is fitted to a vertex. All pairs that 
pass a cut on the vertex-fit mass {AitIt^+t^- = \m^+^- — mxo\ < 10.8 MeV) and a very loose vertex- 
fit probability cut (P(x^) > 10^^) are retained for further consideration. Several more cuts are 
applied to the candidates to reject combinatorial background and mesons that are not B 
decay products. We require the transverse decay-length rjjEC = \J {xk° ~ ^bsY + {Vk^ ~ VbsY be 
between 0.2 and 40.0 cm. Here refers to the fitted vertex position and BS refers to the position 
of the beamspot (the center of the luminous region), which is determined in BABML approximately 
every ten minutes. The inner r^EC requirement removes random vr+vr^ combinations that most 
likely come from a common point (the event primary vertex or a short-lived secondary decay) and 
just happen to have a mass inside the allowed Am^+^- range. Candidates that fail the outer 
^DEC criteria tend to be from calorimeter splash-back tracks that once again just happen to pass 
the Am^+^- cut. We also require that the angle between the transverse flight vector (of which 
fDEC is the magnitude) and the transverse momentum vector of the be less than 200 mrad. 
These cuts, along with the continuum rejection criteria described below, have been optimized in a 
signal-blind study to produce the largest signiflcance for a branching fraction measurement. The 
cuts were optimized using signal MC to represent the B^ —>■ K^K^K^ decay, and on-resonance data 
in mass sidebands, designed to exclude all signal candidates, to model backgrounds in the signal 
m(i^°A'°i^°) region. 

3.2 RECONSTRUCTION 

All combinations of three — > vt+tt^ candidates, where none of the candidates share a charged 
track, are used to form B^ — > K^K^K^ candidates. The momentum is calculated in the 
vr+vr" vertex fit, but because m^+^- is not constrained to rrij^o in the vertex fit, we use Ej^o = 

^JpI+^- + rn\o in place of the energy. 

We use two kinematic variables to separate the B^ K^K^K^ signal from backgrounds. The 
energy difference AE = Eb — \/s/2 is reconstructed from the energy of the B candidate in the e~^e~ 
CM frame Eb and the total energy ^/s. The AE mean value is expected to be near zero for signal 
events, and the AE resolution for signal events is about 18 MeV. The beam-energy-substituted 
mass is defined by rn-Es = y (-5/2 + Pi ■ 'PB)'^/Ef — where {pi,Ei) is the four-momentum of the 
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initial-state e'^e~ system and ps is the momentum of the B candidate, both measured in the 
laboratory frame. The rriEs resolution for signal events is about 2.6 MeV. We retain candidates 
with \AE\ < 0.30 GeV and 5.22 < tties < 5.30 GeV (this is referred to as the bounded AE-rriEs 
plane) . 

Studies of off-resonance data and MC samples have shown that the largest background source 
that passes the above cuts comes from continuum e^e~ uu/dd/ ss/cc events. For this reason we 
choose to cut on three variables that are commonly used to reject continuum background. The most 
powerful continuum-rejection cut is | cos6'j'| < 0.8, where cos Ot is the cosine of the angle between 
the B candidate thrust axis and the thrust axis of the remaining charged tracks and photons in 
the event. The | cos distribution is fairly uniform for signal events, and is peaked near 1 for 
continuum events. We also require —5.0 < J- < -1-1.0, where T \s a, Fisher discriminant [7j based 
on zeroth and second momentum-weighted Legendre polynomial sums of the remaining tracks and 
photons (JT = 0.5264 - 0.1882£o + 0.9417^2 ), and R2 < 0.3, where R2 is the ratio of second to 
zeroth Fox- Wolfram moments j8|. 

After all the above cuts are applied, there are a small (< 1%) number of events with more than 
one 3K^ candidate per event. While the MC appears to properly model this effect, we choose to 
simplify the analysis by using only the candidate with the smallest Z'(Am^+^__)^ in each event, 
where the sum is over the three Kg candidates making up the B^ — > K^K^K^ candidate. 



4 MAXIMUM LIKELIHOOD FIT FOR YIELD 

We take all the combinations that pass the above cuts and extract the yield of signal B^ K^K^K^ 
events (Ns) with an unbinned extended maximum likelihood (ML) fit, where the likelihood is: 



£ = exp (^-E^Jn 



(1) 



the sum over i corresponds to four categories of signal and background, the product over j cor- 
responds to the 508 B° K^K^K^ candidates that pass all the requirements in the previous 
section, and Vij is the probability density function (PDF) for the ith category evaluated for the j'th 
candidate. The PDFs used for the four components of the fit are described in the next sections. 
For all except one category, the PDFs are the products of one-dimensional PDFs in tties and AE: 

Vi = Vi{mEs)Vi{AE). 

The four categories whose yields Ni are determined by the fit are the above-mentioned number of 
signal events Ng, the number of continuum BG events NcbG: the number of BB events in which the 
candidate comes from random combinations of tracks that may or many not be true — > tt^tt~ 
decays (the "non-peaking" BB BG) Nbnop, and events where the three Kg candidates come 
from the same B decay but that B decay is not a signal B^ — > K^K^K^ decay (the "peaking" 
BG) Nbpbg- The default fit requires all Ni > 0, but we loosen this requirement as a systematic 
cross-check, and see that the signal yield changes by a small amount. 

The dominant source of peaking BG in the part of the AE-m-^^ plane over which we do our 
fit is S — > 3-fC^7r decays, where the vr^ or vr'^ is missed. Many of the B branching fractions for 
decays that contribute to this background are known poorly or not at all including the decays 
B^ 2K^K*° and B+ 2K^K*+, which should dominate the peaking BG in the bounded AE- 
rriEs plane. Like the signal B^ — > K^K^K^ events, we remove the peaking BG from the generic 
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BB MC before we fit it to extract the parameters used to describe the non-peaking BB BG PDF. 
We generate large MC samples with only BB peaking BG, extract the peaking PDF from these 
samples, and allow the BB peaking and non-peaking yields to float separately. One reason the \AE\ 
cut is so wide (< 300 MeV) on the bounded AE-uies plane is to allow a large enough region to fit 
the peaking BG. We do not attempt to extract a branching fraction for the peaking BB BG; our 
goal is to determine the size and shape of the peaking BG, so that its presence at an unexpectedly 
large level does not distort the shape of the non-peaking BB BG and lead to a systematic bias on 
Ns- It turns out that the fit to the data requires very little peaking or non-peaking BB BG, and 
the signal yield is changed negligibly if either or both of these populations are fixed at zero. 

4. 1 SIGNAL EFFICIENCY AND PROBABILITY DENSITY FUNCTION FOR 
SIGNAL 

The signal PDF parameters are determined from ML fits to reconstructed and selected signal MC 
events from a sample of 148k generated B^ — > K^K^K^ decays. The signal distribution in AE is 
well described by a double- Gaussian (one for the "core," another for the broader "tail"). The five 
parameters that describe this double- Gaussian PDF are shown in Tabled The signal PDF in niES 
is well described by a bifurcated Gaussian. The values of the three parameters determined from 
the fit to the signal MC are also shown in Tabled Histograms of the signal MC, with the derived 
PDFs shown as overlaid curves, are shown in Fig. [21 A fit to the signal MC with these parameters 
gives a signal yield of Ns = 8147 it 91 events, or a signal efficiency of (5.50 it 0.06)%. 




mEs(GeV) AE (GeV) 



Figure 2: Distributions of mES and AE for selected candidates from MC simulated B Kg Kg Kg 
events. The curves overlaid on the MC data points are the signal PDFs determined from a unbinned 
fit. The parameters from the fit are shown in Table 1. 

4.2 PROBABILITY DENSITY FUNCTION FOR CONTINUUM 

Various studies suggest that a linear function is sufficient to describe the continuum BG in AE, so 
we use that for the AE PDF. For niES the standard PDF is the Argus function 9 , and this proves 



11 



Table 1: Parameters used in the default ML fit. The parameters and their statistical errors 
are derived from MC studies described in the text. The signal PDF is proportional to 
/coreexp( ~^^^^^"°''°^ ) + (1 - fcore) exp( ~^^fj^''"'^ ) , the signal ruEs PDF is proportional to 

core ^ ^"^taii ^ 

6'(/i — m^g)exp( ^^^i ^'^ ) + 0{mEs — A*) exp(— the continuum and non-peaking BB 

left right 

ruES PDFs are proportional to niEs^l - (^^^exp(-e[l - {^?])- 



n if rTiiTTnonPTi"!' 


X CXL CXlliC uCl 


\/ a 1 1 1 p 




ATP n 


D.y It u.z ivie V 




Ai? (Tcore 


14.3 ± 0.2 MeV 






2 ± 1 MeV 


Signal 




37 ± 1 MeV 




Ai? fcore 


(85 ± 1)% 




niES M 


5279.8 ± 0.1 MeV 






2.8 ±0.1 MeV 




(bright 


2.3 ±0.1 MeV 


Common BG 


rriES mo 


5289.8 MeV 


Continuum 


rriES C 


-17 ±8 


BG 


AE linear slope 


-1.9 ±0.3 


Non-peaking 


rriES C 


-47 ± 24 


BB BG 


AE linear slope 


-3.2 ±0.6 



sufficient. There is one parameter (the slope) for the AE PDF and two parameters (mo and ^) for 
the Argus function. The parameters mo and ^ are correlated, and if the mo parameter is set too 
low, there can be problems with the fits. To avoid this problem, mo is fixed to 5.2898 GeV in the 
fit. It is, however, allowed to fioat when systematic errors are evaluated. Both the AE slope and 
^ are fioated in a fit to a continuum MC sample corresponding to 77 fb~^ of integrated luminosity. 
The projections of the MC events and the continuum BG PDFs are shown in Fig.|31 The values of 
the parameters used in these PDFs and in subsequent fits are shown in Tabled 

4.3 PROBABILITY DENSITY FUNCTION FOR NON-PEAKING BB BACK- 
GROUND 

The non-peaking BB BG has a similar source as the continuum BG (random combinations of 
real and fake mesons) and so is expected to have a similar functional form. We used the same 
description (linear function for AE and Argus function for rriEs) as used for the continuum BG and 
fit them to BB MC samples with any candidates identified as signal or peaking BB BG removed. 
We used a mo parameter in common with the continuum BG, but the AE slope and Argus ^ 
parameter are allowed to fioat to different values from those used for continuum BG. The values 
determined by the fit to the non-peaking BB MC samples are also shown in Tabled 
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Figure 3: Distributions of niEs and AE for selected candidates from uu/dd/ss/cc MC simulated 
samples. The curves overlaid on the MC data points are the PDFs used to describe the continuum 
BG in the bounded AE — uies plane. 



4.4 PROBABILITY DENSITY FUNCTION FOR PEAKING BB BACKGROUND 

The distributions for the peaking BB BGs have a very different functional form from those for the 
continuum and non-peaking BB BGs. We have parametrized the peaking BB BG as was done in 
Ref. ^0], by taking the two-dimensional (2D) histogram of all candidates (which pass the analysis 
cuts) from specially generated B^ 2K^K*^ and B~^ — > 2K^K*~^ MC samples and using this 2D 
histogram as the PDF. 

The 2D histogram PDF, with the same binning used in the default fit, is shown in Fig. ^ By 
using a 2D histogram for the PDF, a large correlation between AE and tues, not present for the 
continuum or non-peaking BB BGs, is properly taken into account. Systematic errors due to the 
binned nature of this PDF are discussed later. 



5 FIT FOR YIELD 

We use the PDFs described in the preceding sections and the parameters derived from fits to MC 
samples and listed in Table ^ for our default fit to the on-resonance data sample. The values for 
the populations of the various components of the fit are 

iNs,NcBG,NBNOP,NBPBG) = (71 ± 9, 428ti, Ot^s, 9 ± 8). 

The fit requires no BB non-peaking BG, but does allow for a small but not significant amount 
of BB peaking BG. The projections of the data and the fits on the ruEs and AE axes are shown 
in Fig. [3 The small contribution of the BB peaking BG is the non-overlap of the solid line (all 
components of the fit) and the dashed line (the continuum BG component of the fit) in the region 
AE < -0.10 GeV in the AE plot. 
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Figure 4: The two-dimensional histogram of AE vs. niES used as the PDF for the peaking BB 
BG. 
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Figure 5: Distributions olniES and AE for selected candidates from the on-resonance data sample. 
The solid curve overlaid on the data corresponds to the sum of all PDFs, with their parameters 
(Table and the signal and background fractions returned by the fit. The dashed curve is the 
contribution from the continuum BG. 



6 SYSTEMATIC STUDIES 
6.1 CUT-AND-COUNT ANALYSIS 

We use a simple "cut-and-count" analysis to cross-check the results of the maximum likelihood fit, 
and to study several sources of systematic uncertainty. For the cut-and-count analysis we define 
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a signal region centered on the expected signal in m^o, and AE. The signal region is defined by 
5.2704 <mEs < 5.2884 GeV and -40 < AE < +40 MeV. We define two AE sidebands with the 
same tties cut but with -300 < AE < -100 MeV and +100 < AE < +300 MeV. The sum of 
the number of entries in the sidebands (57) scaled by the ratio of the area in the signal region to 
the area in the sidebands (0.2) gives an estimate of the number of BG events in the signal region 
(11 ± 2), where 78 events are observed, for a signal yield of 67 ± 9 events. The signal efficiency is 
slightly different between the cut-and-count analysis and the ML fit. The MC efficiency-corrected 
yield is 1295^^53 events for the ML fit and 1258 ± 169 events for the cut-and-count analysis. Given 
the different methods for estimating signal and background in the ML fit and the cut-and-count 
analysis, they are in reasonable agreement. 

6.2 SIGNAL EFFICIENCY VARIATION ACROSS DALITZ PLOT 

The cut-and-count analysis allows one to take all the entries in the signal region and plot the m2K0 
distributions for these candidates, and compare them to those predicted by the reconstructed signal 
MC, which was generated with the assumption of non-resonance phase-space (uniform population 
of the Dalitz plot at generation) . The distributions of the reconstructed m2K0 masses are consistent 
with (reconstructed) three-body phase space, but with such a small number of events in the data, 
we cannot rule out resonance production at the level of a few events per resonance in our data 
sample, or other small deviations from phase-space. 
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Figure 6: The folded Dalitz plot for the 78 candidates (black points) that pass the cut-and-count 
analysis and end up in the signal region. The signal/background for this selection is ~ 5.8. The 
yellow-shaded area is the physically-allowed region. 



This is important because the efficiency calculated in different regions of the Dalitz plot is not 
uniform. It varies by more than a factor of two, mostly due to low efficiency for reconstructing 
Kg near some edges of the Dalitz plot. If the parent distribution for the data is the same as the 
MC parent distribution, this is not a problem, but since we do not know this, we need to assign a 
systematic error for the possibility that it is not. 
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We do this by dividing the "folded" Dalitz plot (one of the unique sextants of the 3Kg Dalitz 
plot) into 21 bins and calculating an efficiency for the cut-and-count analysis for each bin individ- 
ually. The folded Dalitz plot we use is shown as the shaded area in Fig. El the points are the 78 
candidates in the signal region of the cut-and-count analysis for the data. This folded Dalitz is 
achieved by ordering the three unique m{2K^) combinations for each K^K^K^ candidate 

niMAxi^K^) > mMiD{2K^) > niMiNi^K'^) and plotting niMiDi^K'^) vs. mMiNi^K^)- 

The entries in the signal region are then individually corrected for efficiency depending on what 
bin they populate, and the yield calculated this way is compared to the yield when all events are 
given the same (average) efficiency. The yield differs by 4.2% between the two ways of calculating 
the efficiency; we take this as the systematic error estimate due to a nonuniform population of the 
Dalitz plot. 

6.3 K° RECONSTRUCTION EFFICIENCY 

There is a small but well-measured disagreement between the — > 7r"'"7r~ reconstruction effi- 
ciency in the data and the one reported by the full detector MC. We correct the efficiency and 
calculate a systematic error on how well we know the reconstruction efficiency. The efficiency 
in the MC simulation and data is measured as a function of transverse decay radius {rDEc)^ 
transverse momentum and polar angle in the BABAR detector, and also for periods with different 
detector running conditions. A correction is calculated for each of the three candidates in a 
reconstructed MC event, and the product of the three correction factors (taken as the B^ candi- 
date correction factor) is averaged over all selected events in the signal MC sample. We do this 
for several different sets of measured efficiencies, each produced with different selection criteria. 
The corrections evaluated for different selection criteria are consistent, and the average correction 
factor is Sdata/^MC = 0.950 ± 0.014. 

The error includes a statistical error for the tables used to calculate the correction, a per- 
charged-track systematic error, and a per-K^ systematic error. The quadrature sum of all these 
error estimates is 10.1%. This is the dominant source of systematic error for this measurement. 

6.4 SIGNAL PARAMETRIZATION 

We allow each of the eight parameters in the signal PDF to float in the fit, one at a time. The 
quadrature sum of the change in signal yield from these eight variations is 4.3%, and we take this as 
the systematic error estimate on our signal parametrization. Since many of the signal parameters 
are correlated, we also perform a fit where six of the signal parameters are free; only the AE 
tail-Gaussian mean and width are fixed to the values derived from the MC. The change in Ns from 
letting all these parameters float together is 2.4%. This is a variant of the above (larger) estimate 
that allows correlations between parameters to be taken into account, but we will use the larger 
estimate as a more conservative estimate. 

6.5 BACKGROUND PARAMETRIZATION 

As with the signal parameters, we allow each of the five BG PDF parameters to float in the flt, 
one at a time. For the binned-histogram peaking BB PDF, we increase and decrease the bin size 
by a factor of two and allow different levels of smoothing of the histogram. We take the largest of 
these variations as the systematic error due the peaking BB PDF, and add it in quadrature with 
the changes in signal yield from letting the BG fit parameters float. The fractional systematic error 
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estimated from all these variations is 0.8%. While the fit will not support all five BG parameters 
floating at once, we take the two parameters with the largest correlation (the A£' slopes for the 
continuum and non-peaking BB BGs) and let them float together. The signal yield changes by 
0.4%, less than the quadrature sum of the two parameters allowed to float separately. While the 
populations of the BG categories change noticeably with all these parameter variations, the total 
BG yield and the signal yield are quite stable. 

6.6 FIT VALIDATION 

We perform studies with parametrized MC simulations in which many samples of the same size (and 
category populations) as the data are generated from the PDFs. We also perform MC studies in 
which the background events are generated from the PDFs but the signal samples were extracted 
from the full-detector MC sample and fit along with the background samples. The means and 
uncertainties for the yields are all consistent with expectations, and no significant corrections for 
biases or systematic error contributions are required. For 2000 fitted toy MC samples, 48% have a 
larger value of — In C than that for the fit to the data. 

There is a 15.9 fb^^ sample of data taken at CM energies just below the T{4S) resonance (the 
off-resonance data sample), which contains no S-meson decays. This data is corrected for a shift 
in the mss endpoint due to different beam energies. The data is subjected to the same selection 
cuts and the same ML fit as the on-resonance data. The signal yield from this fit is consistent with 
zero events. 

6.7 OTHER FIT VARIATIONS 

As a measure of the sensitivity to background parameterizations, we remove the three BG categories 
one at a time in our default ML fit. With cither (or both) of the BB BG PDFs removed, the yield 
changes by very little (< 0.3%). With the continuum BG PDF removed, the signal yield changed 
noticeably (—5.1%), but the likelihood of the no-continuum-BG fit is much worse than that for the 
default fit. That is, neither of the BB BGs (or their combination) does a particularly good job of 
describing the continuum BG, which dominates the bounded A£'-mES plane away from the signal 
region. 

We remove the restriction that each of the four yields be greater than zero and refit. The signal 
yield changes by -|-1.3%, and the likelihood of this fit is only slightly better. 

We add (separately) quadratic terms to the continuum and non-peaking BB AE PDFs and let 
them fioat in the fit; the signal yield changes by —2.2% and —0.6%. We include these variations 
in the systematic uncertainty along with the other systematic errors estimated from floating the 
background parameters above. 

6.8 CANDIDATE SELECTION CRITERIA 

While the candidate selection cuts listed in Section 3 are quite standard and wc expect the MC to 
reproduce them, we estimate a systematic error on each one to account for the fact that the MC 
might not exactly reproduce the data. Where possible, each cut is removed in turn and the change 
in yield for the data is compared with that for the MC. For the few cuts that are significantly 
correlated, both cuts are removed at the same time. The cuts are also tightened by reasonable 
amounts and the change in yield in the data and MC are compared. Various other studies are 
performed on control channels such as B^ D*~Tr~^, D*~ — > -k' , — > K^tt~^'k~ , which has a 
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similar topology and a large enough branching fraction so that it can be reconstructed with minimal 
cuts. The quadrature sum of the systematic error estimated from a variation for each cut is 5.0%, 
where the dominant contributions are from the Am^+^- cut (3.0%), the cosOt cut (2.5%), and the 
i?2 cut (2.1%). 

6.9 SYSTEMATIC ERROR SUMMARY 

There are two other small systematic errors shown in Table |21 not discussed above: the statistical 
error on the MC used to derive the signal efficiency estimate, and the error on the total number of 
BB pairs in our data sample. With these added in quadrature with the systematic error estimates 
described above, the total (fractional) systematic error is 13.1%. 



Table 2: Summary of fractional systematic uncertainties. 



Source 


Estimated from 


Percent Error 


K'^ efficiency 


detector studies 


10.1% 


BG Parametrization 


vary in fit 


2.4% 


Signal Parametrization 


vary in fit 


4.3% 


Candidate Selection Cuts 


cut variations, studies 


5.0% 


Efficiency variation 


cut-and-count analysis 


4.2% 


Signal efficiency 


MC statistics 


1.3% 


BB counting 




1.1% 


Total 




13.1% 



7 PHYSICS RESULTS 

The branching fraction is calculated from the relationship B = Ns/{£N^-^), where Ns = 71 it 9 
is the signal yield from the fit, and e is the product of emc = 5.50% derived from the signal MC 
and Edata/^MC = 95.0%, derived from the efficiency studies. These combine for an efficiency- 
corrected signal yield of 1363^'^Jgy. The data set corresponds to 211 x 10® B^ and 5° decays, and 
we calculate B = (6.5 it 0.8) x 10~® (statistical error only). We assume that the rate for B^B^ and 
B~^B~ production in T(45) decays is equal. The errors on all quantities except for the signal yield 
are included in the systematic error estimate. If we fix the signal yield to zero in our ML fit, the 
difference between the — In £ of this fit and the default fit gives a statistical significance for our 
observation of 15.6(7. 

The sum of the systematic error estimates is given in Table |2I This results in a measurement of 
B = (6.5 ± 0.8 ± 0.8) X 10-6, where the second error is the systematic error estimate. One sigma of 
the total systematic error estimate (not just the ones that pertain to the signal yield) corresponds 
to 9.3 (efficiency-uncorrected) events. When the signal yield is fixed to 9.3 and the data is refit, 
the change in — InC from the default fit corresponds to a significance of 10.9(7. 

We note that this measurement is consistent with the previous Belle measurement, but by itself 
it is more than 2a above the prediction made using B{B^ K~ K^) and the assumption of 

penguin dominance va. B ^ KKK. However, if this difference is confirmed with more data, it may 
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just be evidence of resonant intermediate states occurring at different rates in 5° K^K^K^ and 
((/)-removed) B+ K+K-K+. 

While we have examined the folded Dalitz plot for the cut-and-count analysis and see nothing 
that looks like a narrow resonance (broad resonances cannot be ruled out given the size of the 
data sample), decays like — > XcO^^ may be present and are clearly not part of the relationship 
between K^K^K^ and B^ K~ . To estimate the amount of this type of decay in 

our sample, we reject 2K^ masses within ±50 MeV (about So" of our resolution) of the XcO ^-nd Xc2 
masses. We apply this rejection to the data and the MC samples in our cut-and-count analysis 
and, while a few entries are removed from the data, proportionally slightly more were removed 
from the non-resonant phase-space generated MC, so the efficiency-corrected yield goes up slightly 
(consistent with no change) when the Xc bands are excluded. On the basis of this, we cannot claim 
we observe any contribution to the B^ — > K^K^K^ signal from B^ XcK^, and so we leave the 
inclusive measurement uncorrected. 

8 SUMMARY 

From a sample of 211 x 10^ BB decays recorded with the BABAR detector, we observe a signal of 
71 it 9 B^ — > KgKgKg decays, and with these measure a branching fraction B{B^ K^K^K^) = 
(6.5 =b 0.8 lb 0.8) X 10~^. This result is preliminary. 
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